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Abstract. We present submillimeter observations of rotational transitions of carbon monoxide from J = 2 — » 1 
up to 7 — > 6 for a sample of Asymptotic Giant Branch stars and red supergiants. It is the first time that the high 
transitions J = 6 — - » 5 and 7 — > 6 are included in such a study. With line radiative transfer calculations, we aim 
to determine the mass-loss history of these stars by fitting the CO line intensities. We find that the observed line 
( intensities of the high transitions, including the J = 4 — > 3 transition, are significantly lower than the predicted 

values. We conclude that the physical structure of the outflow of Asymptotic Giant Branch stars is more complex 
than previously thought. In order to understand the observed line intensities and profiles, a physical structure 
with a variable mass-loss rate and/or a gradient in stochastic gas velocity is required. A case study of the AGB 
star WX Psc is performed. We find that the CO line strengths may be explained by variations in mass-loss on 
time scales similar to those observed in the separated arc-like structures observed around post-AGB stars. In 
addition, a gradient in the stochastic velocity may play a role. Until this has been sorted out fully, any mass loss 



Or 



determinations based upon single CO lines will remain suspect. 
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1. Introduction a dust driven wind will develop. The mass-loss rates in- 
jJJ ■ crease from M « 10~ 7 to a few times 10~ 5 M Q yr _1 , while 
Low and intermediate mass stars (1 < M < 8M ) end the AGB star evolves from th e Mira phase to an OH/IR 
their life on the red pant branch and asymptotic giant gtar || van d er Veen fc HabindllQsl) . Recently, it has been 
branch (AGB; see lHabing] | l996l and references herein). suggested that higher masS -lo SS rates can be achieved for 
During the AGB phase, the stars have very extended ten- oxygen _ rich AGB sta rs. I.Iusttanont et all [l996) find that 
uous atmospheres and shed almost their entire hydrogen- 0H 26 .5+0.6 has undergone a recent increase in mass loss, 
rich envelope through a dense and dusty stellar wind. In leading to ft curren t rate of5.5- lQ- 4 M yr" 1 , a result re- 
case of OH/IR stars, mass-loss rates can be so high that C£mtly confirmed hv lFong et a,lJ {zool . Even higher mass- 
the dust shell completely obscures the central star, and logs rateg were found for another oxyge n-rich AGB star, 
the object is observable only at infrared wavelengths and IRAS 1634 2_3814, for wh ich the mass-lo ss rate may be as 
through molecular line emission at radio wavelengths. The high ag „ 1Q - 3 Mq yr -i fo iikstra et alJliffli . A similar 
AGB phase is one of the few occasions m stellar evolution rate of a few timeg 1Q - 3 Mq -i ig found fo r the carbon _ 
when time scales are not driven by nuclear (shell) burn- rich evolyed gtar AFGL 268g jskinner et alJliS . 
ing but by surface mass loss. Helped by the low surface 

gravity and strong stellar pulsations, gas can move away AGB stars are important contributors of dust to the 

from the star and will gradually cool. When the temper- interstellar medium (ISM); it is estimated that a substan- 

ature drops below ~1400 K, dust formation occurs, and tial fraction of the i nterstellar du st is produced by oxygen- 

rich AGB stars (e.g. lGehrJl989|) . In the outflow of evolved 

Send offprint requests to: F. Kemper (kemper@astro.ucla.edu) stars with an oxygen-rich chemistry the dust composition 

* SIRTF Fellow is dominated by silicates, both amorphous and crystalline 
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(e.g. ISvlvester et~aHll999t iMolster et afl 12002^ . The ap- 



pearance of crystalline silicate features in the far-infrared 
spectra of AGB stars seems to be correlated with a high 
optical depth in the amorphous silic ate resonance at 9.7 



Aim and hence a high mass-loss rate ( 


Waters et al.lll996l 


ICami et al.lll998HSvlvester et al.lll99£ 


). This could be in- 



terpreted as evidence that a certain threshold value for the 
density is required t o form crystalline silicates. However, 
iKemner et al showed that observational selection 

effects may play an important role in detecting crystalline 
silicates in AGB stars with low mass-loss rates. Therefore, 
the relation between mass-loss rate and crystallinity re- 
mains unclear at present. 

In order to further study the correlation between the 
wind density and the dust composition, reliable mass- 
loss rates should be determined. Mass-loss rates of AGB 
stars can be obtained from the thermal emission from 
dust , predominan tly coming from the warm inner regions 
(e.g. lBediinlll987IK They can also be inferred from obser- 
vatio ns of molecular transit ions, in particular from CO 
/e.g. iKnann fc Morris! 1 198.^ . A catalogue compiled by 
ILoup et al.1 (jl993^> lists observations of the CO J = 1 — > 
and J = 2 — > 1 transitions of both O-rich and C-rich 
AGB stars. (Hereafter we will use for these rotational tran- 
sitions the notation CO(l— 0) etc.) The mass-loss rates 
of a large number of objects from the catalogue are de- 
rived. However, the derived mass-loss rates seem to be 
underesti mated for OH/IR st ars, compared to the dust 
mass loss. iHeske et al.l l)l990|) have studied the correla- 
tion between IRAS colours and mass-loss rates derived 
from CO(2— 1) and CO(l— 0) observations. In the case 
of very massive dust shells, they find that the intensity 
of the CO(l— 0) transition is too low compared to the 
CO(2— 1) transition, which they suspect to be due to a 
mass-loss rate increase over time. This then hints towards 
a superwind phase, which is generally believed to be im- 
portant in the evolution of a M ira towards an OH/IR 
star (e.g. llben Jr. fc Renzinilll983L and references herein). 
The superwind model was initially introduced to explain 
the amount of mass seen in planetary nebulae a ssuming 
that Miras are the progenitors of these nebulae i|Renzinil 
Miras are believed to evolve into OH/IR stars when 
they suddenly increase their mass-loss rate with a factor 
of -100. 

As the inner regions are warmer they are better probed 
by higher rotational transitions. Thus a sudden density 
jump sh ould be detectable in the CO lines. Model calcula- 
tions bv l.lusttanont et al.l l|l99ft ) have demonstrated this 
effect for OH 26.5+0.6, using observations of rotational 
transitions up to CO (4— 3). Unfortunately this transition 
is not sufficiently high to firmly establish the recent on- 
set of a superwind, as its excitation t e mper ature is only 
55 K. Nevertheless, l.lusttanont et al.1 l|l996h found that 
the peak intensities of these lines were significantly higher 
than what could be expected based on the extrapolation 
of the observed line strength of the CO (2—1) transition 
and the upper limit obtained for the CO(l— 0) transition, 
assuming a constant mass-loss rate. Similar results are 



Table 1. Technical details of the JCMT heterodyne re- 
ceivers. The columns list the used receivers, the frequency 
windows at which they operate, the observable CO ro- 
tational transition, the beam efficiency ?y m b and the half 
power beam width (HPBW). 



receiver 


Frequency 


CO transition 


7?mb 


HPBW 




(GHz) 








A3 


215-275 


CO(2-l) 


0.69 


19.7" 


B3 


315-373 


00(3-2) 


0.63 


13.2" 


W/C 


430-510 


00(4-3) 


0.52 


10.8" 


W/D 


630-710 


00(6-5) 


0.30 


8.0" 


E 


790-840 


00(7-6) 


0.24 


6.0" 



reported for other AGB stars (e.g. iGroenewegenl Il994lj: 
iDelfosse et al1ll997^ . 

The work presented here aims to determine the mass- 
loss history of a number of oxygen-rich AGB stars with 
an intermediate or high optical depth in the near- and 
mid-infrared. For the first time, observations of rotational 
transitions up to CO (7—6) have been obtained (T cx = 
155 K) which probe the more recent mass-loss phases. In 
Sect. [5] we describe the observations and data analysis. 
Sect. |21 describes the model. Our results are discussed in 
Sect.0] Concluding remarks and an outlook to future work 
is presented in Sect. [3] 

2. Observations and data reduction 

2,1. Instrumental set-up 

Observations of the 12 CO(2-l), (3-2), (4-3), (6-5) and 
(7—6) rotational transitions in the outflow of evolved 
stars were obtained during several observing periods be- 
tween April 2000 and September 2002 using the James 
Clerk Maxwell Telescope (JCMT) on Mauna Kea, Hawaii. 
For this purpose, all five different heterodyne receivers 
available at the JCMT were used, including the new 
MPIfR/SRON E-band receiver which operates in the 
790-840 GHz frequency range. A description of this new 
receiver is given in Sect. 12.21 The technical details and 
beam properties of the JCMT set up with the appro- 
priate heterodyne receivers are summarized in Table Q 
Observations with the B3- and W-receivers were per- 
formed in double sideband (DSB) and dual polarization 
mode. The DSB mode was also used for the observations 
with the MPIfR/SRON E-band receiver. The bandwidth 
configuration of the receiver, and hence the spectral reso- 
lution was determined by the expected line width of the 
CO lines. We used bandwidths of at least twice the ex- 
pected line width to have a sufficiently broad region for 
baseline subtraction. Estimates for the line width - which 
is determined by the outflow velocity - were based on pub- 
lished values of li ne widths of the CO(l— 0) transition (e.g. 
ILoup et aflll993L and references herein). 

We used the beam-switching technique to eliminate 
the background. The secondary mirror was chopped in 
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azimuthal direction over an angle of 120". Over these small 
angles the noise from the sky is assumed to be constant. In 
case of extended sources we used a beam-switch of 180". 



2.2. The MPIfR/SRON 800 GHz receiver 

The observations of the CO(7— 6) line were made with the 
MPIfR/SRON 800 GHz receiver in October 2001. This 
PI system is in operation at the JCMT Cassegrain fo- 
cus cabin since spring 2000. The receiver consists of a 
single-channel fixed-tuned waveguide mixer with a diag- 
onal horn. The mixer consists of a Nb SIS junction with 
NbTiN and Al wiring layers fabricated at the University 
of Groningen, The Netherlands. De tails on the fab r icatio n 
of similar devices can be found in Ijackson et alJ l|2000|) . 
Measured receiver temperatures at the cryostat window 
are Tr^ ~ 550 K DSB. The receiver has an intermediate 
frequency of 2.5 — 4 GHz. System temperatures including 
atmospheric losses varied between 6000-14000 K (SSB) 
at the time of the observations. The beam shape and effi- 
ciency have been determined through observations of Mars 
and yield a deconvolved half power beam width (HPBW) 
of 6" and a main beam efficiency f? m b of 24%. 

2.3. Observations and data reduction 

Our sample of evolved stars is given in Table EU which 
also indicates the distances towards the programme stars. 
The sample includes AGB stars and red supergiants. In 
Table lA~2l an overview of the observed transitions is given, 
including cumulative integration times and the observing 
date. The data were obtained over a long period from April 
2000 until September 2002 in flexible observing mode, and 
are part of a larger ongoing programme. During the ob- 
servations, spectra of CO spectral standards used at the 
JCMT were also obtained. If necessary, a multiplication 
factor was applied to the observations of our sample stars, 
to correct for variations in the atmospheric conditions. 
These factors are listed in Col. 4 of Table I A. 21 and are 
based on measured standard spectra. Reliable standards 
are only available for the transitions observed with the 
A3-, B3- and W/C-receivers, for which the flux calibration 
accuracy is around 10%. For the W/D- and MPIfR/SRON 
E-band reliable standards for our lines of interest are lack- 
ing. Therefore we estimate that the absolute flux calibra- 
tion in these bands has an accuracy of 30%. 

Tablcn]lists the beam efficiencies ?y m b for all receivers. 
The main beam temperatures were calculated according 
to T m b = Tj^/rimb) where T£ is the measured antenna 
temperature. These main beam temperatures can directly 
be compared to observations from other telescopes. 

The reduced data is presented in Table lA~3l A linear 
baseline has been subtracted from the raw data, and the 
spectrum has been rebinned to improve the signal-to-noise 
ratio. We aimed to cover the line profile with at least ~80 
bins, which limits the rebinning factor. The bin sizes after 
rcbinning and the corresponding r.m.s. values arc listed in 
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Fig. 1. Correction of the profile of the CO(3— 2) transi- 
tion of VX Sgr. The dotted line represents the observa- 
tion in which the interstellar contribution is clearly vis- 
ible. Ignoring the interstellar contribution results in the 
solid line, which is used to obtain the integrated intensity. 



Cols. 4 and 5 of Table lA~3l Emission lines were detected 
in almost all observations, except for a Sco CO (3— 2) and 
OH 104.9+2.4 CO(6-5) and (7-6), for which we only 
obtained upper limits on the main beam temperatures. 
The line profiles of all transitions are shown in Figs. IA.il - 
IA.21I In some cases, interstellar lines are visible in the 
spectrum, for example in VX Sgr. To determine the in- 
tegrated intensities we have cut the interstellar lines out 
of the spectrum, and interpolated both parts of the spec- 
trum, as is demonstrated in Fig. ^ The resulting profile 
was integrated to obtain I, which is the integrated inten- 
sity in K km s . The system velocity Vlsr & n d the ter- 
minal expansion velocity Voo are estimated directly from 
the line profile. The lines show a wide variety of shapes. 
There are parabolic line profiles, like those of WX Psc 
fFig. rOJ . IRC+ 50137 (Fig.EOll. AFGL 5379 (Fig.jOJjl, 
CRL 2199 (Fig. lA~T2)l and OH 104.9+2.4 (Fig. lA~2l) . 
These paraboli c line profiles indicate that the lines are 
optically thick l|Morrislfl980t) . On the other hand, many 
objects show signs the double-horned profiles indicative 
of an optically thin molecular layer. The most illustrative 
example is VY CMa, in which the CO(2-l) and CO(3-2) 
transitions clearly show a double-peaked profile, although 
the peak around the central velocity indicates a more com- 
plex outflow structure (Fig. IA.9")l . In addition, some flat- 
topped profiles are observed, most notably those of VX Sgr 
(Fig. IA.11|) . These flat-topped line profiles are considered 
to be characteristic of molecular l ayers which have r ~ 1 
at these frequencies (MomsJ li980h . 

Ijusttanont et all 1119961) have observed OH 26.5+0.6 
with the JCMT as well and report that they find line 
intensities I = 25.8 and 36.0 K km s" 1 for the CO(3-2) 
and (4—3) transition respectively. In addition they have 
scaled IRAM observations of the CO(2— 1) to a 15m dish, 
to mimic the JCMT. The intensity of this line turned out 
to be 7.8 K km s _1 . Their results agree well with our 
results in case of the CO(2— 1) and (3—2) transition, but 
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they have observed an intensity of a factor of ~2 higher 
for the CO (4— 3) transition. The origin of the discrepancy 
with our results is unknown. 

3. Physical conditions in the outflow: a model 

The observed line profiles provide information on the 
physical structure of the outflow of these AGB stars, 
as the spectral resolution at the observed frequencies is 
sufficiently high to resolve the velocity structure. The 
terminal expansion velocity can be derived directly 
from the width of the line profile (Table IA.3(1 . The 
model w e use to analyze the CO data is based on a 
study by Schonberd lll988(l and was previously used by 
I.Tusttanont et al.l l|l994l) ."The interpretation of our obser- 
vations using this model is discussed in Sect.rj] 



3.1. Description of the model 

The code consists of two parts: The first part solves 
the radiation transfe r equation in the co-moving frame 
(Mi halas et alJ Il975h . computes the level populations 
(in full non-LTE) and iterates until level populations 
and radiation field are consistent. For solving the 
level populations, a Newton-Raphson method is used 
dSchonberg fc Hemne! Il98ri . The calculations take into 
account (de-) excitation through collisions, of which the 
rate is defined by the thermal velocity distribution, calcu- 
lated from the local temperature, as well as (de-)excitation 
induced by a local radiation field and spontaneous de- 
excitation. The code treats pure rotational transitions in 
the ground and first vibrational levels, which are con- 
nected through these collisional and radiative transitions. 
The model can calculate the populations of as many as 
50 levels at once, and is also applicable to molecules other 
than CO. The non-LTE rate equations to determine the 
level populations are described by 



J2n j (A ji + C ji + B ji J ij ) = 



(1) 



A change from level i to level j can be induced by 
collisional transitions (with the collisional rate Cy ) and 
radiative transitions, including spontaneous emission (^4y, 
where Ajj = for i < j) and stimulated emission and 
absorption (BijJij). The collisional transition rates Cjj 
are taken fr om laboratory measurem ents and potentials 
calculations l|Flower &: Lau nav 1985) and are extended up 
to J = 30. 

The line profile integrated mean intensity Jy- consists 
of two components: 

i) The continuum radiation, originating from dust locally 
present. This radiation field can be switched off, by 
assuming there is no dust present in the considered 
part of the outflow. 



ii) Line radiation originating from a local region. The 
size of th is regio n is defin ed by a velocity which 
ISchonberd (|l988h and also I.Tusttanont et alJ l|l994h 
have referred to as stochastic velocity v s to- The nature 
of this stochastic velocity is not specified, but physi- 
cally should consist of a thermal component i>therm and 
a turbulent component fturb, given by 



Vsto 



Vth, 



+ (Vtu 



») 2 



(2) 



In the outflow, the stochastic velocity is assumed to be 
constant and in almost all cases dominated by turbu- 
lence. The effect of the stochastic velocity is Doppler 
broadening of the lines, which is taken into account in 
the radiative transfer. 

In the second part of the code, the calculated level 
populations are used as input to determine the observable 
line profiles by ray-tracing. Again the stochastic velocity 
is used, this time to determine the width of the interaction 
region along the line-of-sight to the observer. Integration 
over the full beam, for which the telescope parameters are 
required, yields the emergent line profile. 



3.2. Free parameters 

The model has a number of free parameters (see Table FJJ). 
In this section we will discuss the various parameters and 
their relevance for the model calculations. 



3.2.1. Density profile 

The density profile p(r) of the outflow determines the col- 
lision probabilities and optical depths required to solve 
Eqs. (TJJ and to calculate the line profiles. The density 
profiles follows from the equation of mass continuity 



p(r) 



M 



(3) 



4-7rr 2 v cxp (r) 

where the expansion (or outflow) velocity profile used in 
the model is defined by 



b 

1 - - 

r 



(4) 



In this equation represents the terminal velocity. 
Constant b is chosen such that the expansion velocity at 
the stellar surface is given by v cxp (i?») = lO -2 ?;^. The 
density structure is set by the following input parameters 

i) The gas mass-loss rate M determines the mass input at 
the inner radius of the circumstellar shell. Our model 
allows us to simulate the effect of a time- variable mass- 
loss rate introducing one jump in the mass-loss history 



at an arbitrary point in the outflow (i 



supcrwind j ? 



where 



the density can increase or decrease with a specified 
factor. Except for this jump the mass-loss rate is con- 
stant, and therefore the density profile scales with the 
current mass-loss rate at r < r supcrw i n d and with the 
past mass- loss rate at r > r supe rwind- 
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ii) The density profile also scales with the outflow velocity 
profile given in Eq. flljl. which is fixed by the terminal 
velocity v^. 

iii) The stellar radius R t determines the base of the wind. 
The density p(i?») at the inner radius follows from 

M and v cxp (i?* ) using the equation of mass continuity 

m 

iv) The outer radius R ou t denotes the extent of the out- 
flow. 

3.2.2. Temperature profile 

The temperature profile T{r) is another important pa- 
rameter that influences the level populations in the cir- 
cumstellar CO, by means of collisions. The tempera- 
ture profile may be compiled self-consistently, i.e. based 
on calculations of realistic heatin g and cooling pro cesses 
/e.g. lOoldreicn fc Scovillel Il97rl j.Tusttanont et al] Il994l 
IChen fc Neufeldl 11993 IZnbko fc ElitzuHl2(WL As a first 
order estimate we have used a power law of the form 
T(r) oc r~ a , where the index a depends on the mass-loss 
rate and is derived from th e outer regions of the te mper- 
ature profiles calculated bv lJusttanont et alJ l|1994|) . 

3.2.3. Dust-to-gas ratio and dust properties 

Unfortunately it is difficult to study the gas and dust 
mass-loss rate completely independent from each other, 
as continuum emission from dust may have an effect on 
the (de-)excitation rates, as described in Sect. 13.11 In par- 
ticular, infrared photons at 4.6 fim pump CO molecules 
from the ground vibratio nal state v = to the first v i- 
brational level v = 1 (e.g. lMorridll98Cl ISchonberdll98S|) . 
The molecules will eventually de-excite to the vibrational 
ground level, but not necessarily to the same rotational 
ground level. This causes a higher population of the higher 
CO rotational levels than which reflects the kinetic tem- 
perature of the gas and the line radiation field. As the 
source of the 4.6 /im radiation is predominantly thermal 
dust emission, the dust-to-gas ratio and the dust opac- 
ity are important input parameters. For simplicity, we as- 
sumed that there was no dust present in the outflows. For 
some of the calculations we did include dust to study the 
effect on the line strengths. In those cases we used a dust 
opacity corresponding to the mixture of solid state com- 
ponents derived for O H 127.8+0.0, a typical OH/IR star 
I Kemper et al.l Eo02j) . The same power-law temperature 
distribution as for the gas is used to calculate the thermal 
emission from the grains, although this is most likely not 
true. 

3.2.4. Velocity field 

The velocity field has already been mentioned as a con- 
straint for the density structure, but it also plays an im- 
portant role in the formation of line profiles. The outflow 
velocity profile (constrained by the terminal velocity Voc 



" 40 



£ 200 



a Sco 


T T Cet 
I 


V438 Oph 


standard 

X 

X 

x x 


fj, Cep 

' \ 


OH104.9 + 2.4 

s 


IRAS21554 


OH44.8-2.3 

f 


OH26.5 + 0.6 10H127.8 + 0.ofCRL 2199 1 IRC +50137 

1 1 1 ^ 

o | i f 

$ i * i t 


a Ori 

<& 


WX Psc 

s§ 1 


AFGL 5379 
§ 


o Cet 


VX Sgr 




IRC +10420 




NML Cyg 

5 

o * j 


VY CMa <» 

V 



►.mo oio> (J oi Giro oioi 

Fig. 2. Overview of integrated intensities for each line ob- 
served in our programme stars. The horizontal axis of 
each panel lists the rotational transitions observed, where 
the spacing between the tick marks is proportional to 
the difference in frequency. On the vertical axis the in- 
tegrated intensity (K km s _1 ) is given. The diamonds 
represent the measured values; in addition the error bars 
are shown (data from Table lA~3l) . Note that only the line 
strengths of VY CMa increase with higher rotational tran- 
sitions. For most of the other stars (except OH 127.8+0.0, 
IRC +50137 and IRAS 21554+6204) CO (3-2) is the 
brightest line. In the upper right corner the relative values 
for the standard model (see Tabled Fig.[3J) are presented 
for comparison (indicated with x symbols). 

and the velocity law given in Eq. Q) and the stochastic 
velocity v sto determine the location and extent of the in- 
teraction regions. As said before, the stochastic velocity is 
assumed to be constant throughout the dust shell. 

3.2.5. Distance and telescope parameters 

The resulting main beam temperatures depend on the dis- 
tance towards the object. In addition, the telescope beam 
size is important to determine what part of the object falls 
inside the beam. In case the circumstellar shell is resolved, 
the pointing displacement (usually 0") should be known 
as well. 

4. Analysis of the results 

Here we will analyze the observations using the model de- 
scribed in Sect. El Fig. shows the intensities integrated 
over line width of each observed line for all our sample 
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Table 2. Parameters of the standard AGB model 
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Fig. 3. Line profiles calculated for a standard AGB model, 
folded with the JCMT beams. The model parameters are 
given in Table [3 



stars. For all sources, except VY CMa, the integrated in- 
tensity increases from the CO(2— 1) to (3—2) transition, 
and decreases again for higher transitions. This is also vis- 
ible in the peak main beam temperatures T m b (Table lA~3)l . 
Since most studies concentrate on the lower transitions 
(up to CO(3— 2)) this was not noticed before. An ex- 
ception ar e the JCMT observation s of OH 26.5+0.6 per- 
formed bv lJusttanont et all l|l996(l . where the CO(4-3) 
transition is included as well. However, as pointed out in 
Sect. 12.31 they observed an increasing line strength with 
increasing rotational transition. This differs from our ob- 
servations of this object that the CO(3— 2) is the strongest 
emission line. 

In order to explain the observational trends, we have 
constructed a standard model assuming physical parame- 
ters widely used for AGB outflows (see Table EJ- We used 
a mass-loss rate of M = 10~ 5 Mq yr _1 , and calculated the 
level populations of the CO gas between the stellar radius 
R* = 4.0 • 10 13 cm and the outer radius R out = 6000 i?*. 



For the terminal velocity we used v a 



15.0 km s" 



and the turbulent velocity was assumed to be w sto = 1.0 
km s . A power-law temperature profile was chosen: 
T(r) = 2000 (r/R„)~ - 7 K. We used for the relative abun- 
dance of the CO gas with respect to molecular hydrogen 
[CO]/[H2] = 3.0 • 10~ 4 , and we ignored the contribution 
of thermal emission from dust to the local radiation field. 
Finally, we placed this system at a distance of 1000 pc, 
and used the JCMT telescope parameters to calculate the 
emerging line profiles (Tabled Fig. 0. The lines show 
increasing peak and integrated intensities with increasing 
line strengths, up to CO (6— 5). The CO (7— 6) line is again 
much weaker which can be explained by the relatively nar- 
row HPBW of the E-band (Table UJ). This transition is 
comparable in strength to the CO(2— 1) transition, for this 
standard set of parameters. This is a general characteristic 
of all other studies calculating the line intensities for corn- 



parameter 



value 



distance 

Woo 
Vsto 
Rin 
Rout 

R, 
M 
T(r) 

[CO]/[H 2 ] 



1.0 kpc 
15.0 km s _1 
1.00 km s _1 
5 R* 
6000 R, 
4.0 • 10 13 cm 
10~ 5 M yr" 1 
2000 (r/R,)-°- 
3.0 • 10~ 4 



K 



dust-to-gas ratio 0% 



monly used AGB para meters (e.g. lGroeneweger]ll994albl : 
lJusttanont etajlliflfll. 

In the following sections we will try to find a set of 
parameters to explain our observations: in general the 
CO(3— 2) is the strongest line, which contradicts the re- 
sults of the standard model. In order to study as many 
stars as possible in a systematic way, we will use a 
line ratio diagram based on the CO(3— 2)/CO(2— 1) and 
CO(4— 3)/CO(2— 1) ratios of integrated intensities, rather 
than trying to fit the intensities and line profiles. 

The two low mass-loss rate AGB stars RV Boo and 
X Her are added to the sample; these stars are the only 
ones for which sufficient reliable line ratios of interest can 
be derived from published JCMT data (see Table IA.4|) . 
For RV Boo the ratios of the integrated intensities are 1.5 
and 2.0 for C O(3-2)/CO(2-l) and CO(4-3)/CO(2-l) 
respectively l(Bergman et al.l l200f)fl . For X Her these 
numbers are 2.1 and 3.4 respect ively ijKnapp et alJ 
Il9 98t iKerschbaum fc Olofssonl Il999h . The observations 



of Knapp et al.1 l|l998|) were obtained with the CalTech 
Submillimeter Observato ry (CSO) and are rescaled to 
the JCMT observations of lKerschbaum fc Olofssonl (^.999) 
such that the line intensities of the CO(3— 2) transitions 
both reflect the same dish size and can be compared to our 
CO(3— 2) and CO(2— 1) observations. One should bear in 
mind however, that RV Boo and X Her ar e not representa- 
tive of AGB stars with spherical outflows. iKahane fc Jural 
(1996) have mapped X Her in CO lines and conclude that 
in addition to a slow spherically expanding shell there 
are indications for bipolar outflows with a higher veloc- 
ity, which carry a signific ant fraction of the ejecte d gas. 
This result is confirmed bv lKerschbaum et alJ l|2003|) . who 
in addition present SiO line observations indicative of a 
circumstellar rotating disk. Interferometric CO line ob- 
servations of RV Boo indicate that this o bject also has a 
disk, possibly showing Keplerian rotation l|Bergman et alJ 
2000|). Therefore, comparison of these stars with our data 
and analysis should be done with some reservation. 

In the literature, we found a sample of six Miras, which 
were observed in all three lines discussed here, using CSO 
(see Table |A~4|) . It is possible to scale these observations 
to the JCMT observations by accounting for the dish size. 
However, we have chosen not to do this, because it is un- 
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Fig. 4. Line ratio diagram. On the horizontal axis the ra- 
tio of the integrated intensities of CO(3— 2)/CO(2— 1) is 
given, while the CO(4— 3) /CO(2— 1) ratio is plotted on the 
vertical axis. The diamonds represent the positions of our 
sample stars, complemented with literature data for X Her 
and RV Boo. Equal ratios are indicated with the dashed 
line. In case of our sample stars, the CO(3— 2) transition is 
stronger than the (4—3) transition, therefore all observa- 
tions can be found in the lower right half of the diagram. 
The only exception is OH 127.8+0.0 which is found in the 
upper left half. The asterisks mark the positions of model 
calculations, where we used the standard parameters (see 
Table |2l . Only the mass- loss rate was varied and is given 
in units of Mq yr _1 . 

known how reliable the rescaled data still is, as little is 
known about the beam filling factor of the various transi- 
tions, while the beam sizes of the telescopes are very differ- 
ent. Instead we chose to compare these observations with 
our model calculations, as will be discussed in Sect. 14. ll 

4.1. A constant mass-loss rate? 

In a line-ratio diagram (Fig. 0}, the observed values oc- 
cupy the lower right half of the diagram, corresponding 
to the region where the CO(3— 2) line is stronger than the 
CO(4-3) line. The values corresponding to OH 127.8+0.0 
are an exception and are found in the upper left half. This 
data point should be treated with care though, as the de- 
tected lines suffer from interference with interstellar ab- 
sorption (Fig. IA.3(1 and therefore the line intensities are 
not well known (see Table IA.3JI . The literature data of 
RV Boo and X Her are also located in the upper left half 
of the diagram. 

First, we assume that the mass- loss rate is constant. 
For five different mass-loss rates (10 -8 , 10~ 7 , 10~ 6 , 10~ 5 
and 10 -4 Mq yr _1 ) we have calculated the emerging line 
profiles, thus cove ring the full r ange in M from Miras 
to O H/IR stars l|Bediinl Il987t Ivan der Veen fc Habind 
1988|. All other input parameters were assumed to have 
the standard values given in Table The predicted 
CO(3-2)/CO(2-l) and CO(4-3)/CO(2-l) ratios of the 




2 4 6 8 10 

co 3-a / co a-i 



Fig. 5. Line ratio diagram for CSO data. On the 
horizontal axis the ratio of the integrated inten- 
sities of CO(3-2)/CO(2-l) is given, while the 
CO(4-3)/CO(2-l) ratio is plotted on the vertical 
axis. The diamonds represent the ratios obtained for a 
sample of six Miras The CO(4-3 ) and CO (3-2) data 
are obtained from lYound (Il995l). and th e CO (2-l) 
data is taken from a study bv iKnapp et alJ l)l998|) . For 
R Hya and \ Cyg additiona l CSO CO(3 — 2) data are 
availa ble, from IKnapp et alJ l|l998l) and IStanek et alJ 
l|l995|) respectively. In the di agram these measu rements 
are indicat ed with 'R Hya K' ijKnapp et al.lll99^) and an 
'X Cvg S' Iptanek et al.lll995j) . Equal ratios are indicated 
with the dashed line. The asterisks mark the positions 
of model calculations, performed for the CSO beam and 
dish size, where we used the standard parameters (see 
Table Only the mass-loss rate was varied and is given 
in units of Mq yr^ 1 . 

integrated intensities were compared to the observed ra- 
tios. 

The model calculations (marked with asterisks) are 
found in the upper left half of the diagram where the 
CO(4— 3) line is stronger than the CO(3— 2) line, and are 
therefore not consistent with the observed line ratios. All 
model line ratios are found in a narrow range to one end of 
the region where the observations are found (see Fig. Q . 
Only the observations of RV Boo match the modelled line 
ratios, but this could be merely a coincidence as RV Boo 
is not a typical AGB star. 

There are not many reports in the literature of AGB 
stars observed in these three lines with the JCMT, but 
we compared the results discussed here with observations 
performed using the CSO. For that purpose, we have re- 
calculated the model line ratios for the CSO beam and 
dish size. A sample of six Miras is consistently observed 
with CSO, where the CO(4-3) and CO(3-2) measure- 
ments are obta ined bv|Youne| lll995 | ) . and the CO(2-l) 
observations bv lKnapp et alJ l|l998). Additional CSO ob- 
servations of t he CO(3— 2) line in two of these objects were 
also included l|Knapp et al.ll998HStanek et al.lll995|) . The 
results are shown in Fig. [S] Similar to Fig. the line ra- 
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tios derived from the standard calculations are found just 
above the dashed line, indicating that the CO (4— 3) line 
should be stronger than the CO (3— 2) transition. However, 
most observations are found well below the dashed line, 
where the CO (3— 2) is the strongest line. \ Cyg falls in 
this region as well. H owever, if we use the measurement of 
IStanek et ail l|l995|) for the CO(3-2) line, the line ratios 
become such that it is found in the same region as the 
model ratios. Possibly this point is unreliable, as it does 
not come from a consistent data set. R Hya seems to be 
an outlier for both CO(3— 2) measurements. Another re- 
markable observation is that the observed CO(2— 1) lines 
seem to be weaker than what is expected from the model 
calculations, given the fact that the calculated ratios are 
closer to the origin of the plot. We may conclude that 
in general the CSO observations occupy more or less the 
same region of the plot with respect to the model ratios as 
our observations. Therefore, in the remainder of this paper 
we will limit our detailed analysis to our JCMT data. 

Apparently, variations of the mass-loss rate alone do 
not change the line ratios enough to significantly increase 
the strength of the CO(3— 2) line with respect to the 
CO (4— 3) line. In the next section, we will investigate to 
what extent variations in the other parameters can shift 
the model calculations such that the line strength ratios 
more closely resemble the observed values. 



4.2. Exploring parameter space 
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Fig. 6. Mosaic of diagrams representing the 
CO(4— 3)/CO(2— 1) ratio on the vertical axis versus 
the CO(3-2)/CO(2-l) ratio on the horizontal axis. We 
used the standard model described in Sect. 0] and Table [21 
and varied for each panel one of the parameters. From 
the upper left corner turning clockwise the investigated 
parameters are: the stochastic velocity i> sto (km s _1 ), 
the outer radius i? ou t the dust-to-gas ratio and the 

distance D (pc). The line ratios resulting from the model 
calculations are marked with asterisks and the observed 
line ratios with diamonds. The dashed lines indicate equal 
line ratios. Note that the ranges plotted on the axes are 
smaller than the ranges in Fig. 0] to improve readability. 



To further explore parameter space, we opted to vary the 
input parameters of the standard model (Table |2J) one by 
one, and compare the line ratios with the observations. 
Combining the changes in line ratios from variations in 
the individual parameters then provides a feeling for the 
range in line ratios that can be covered, and may show 
whether or not it is possible to explain the observed line 
ratios at all. Of course, once a satisfactory match in line 
ratios is achieved by combining the effects of changes in 
individual parameters, fine tuning should be performed to 
fit the observed data in detail. This is necessary as some 
of these parameters might not be completely independent 
from each other, and the precise combined effect on the 
line profile is difficult to predict. 

In Fig. a mosaic of line-ratio diagrams is shown, in 
which the effects of changes in the stochastic velocity, 
the outer radius, the dust-to-gas ratio and the distance 
are shown. In general the effects due to changes in these 
parameters are small. To keep the plots readable, only 
small parts of the original line-ratio diagram (Fig. are 
shown. The modelled line ratios for which these parame- 
ters are varied scatter mainly closely around the observed 
values for RV Boo. In all these modelled line strengths, 
the CO(3-2) line is still weaker than the CO(4-3) line. 
Of course varying the parameters mentioned here causes 
changes in the absolute line strengths, but the line ratios 
are not so much affected. 



In the models where the distance was varied, we placed 
the object progressively closer to the observer, such that 
the beam filling factor is initially less than unity, but in- 
creases with decreasing distance. Although the beam size 
corresponding to the CO(2— 1) transition is larger than 
that corresponding to the CO(3— 2) transition, the line 
formation region of the CO(2— 1) transition is located 
so much further out that the object is first resolved for 
the CO(2— 1) transition. This implies that less emission 
from this line is received by the telescope. When this hap- 
pens, the line ratios increase. The CO(4— 3)/CO(2— 1) line 
ratio increases faster for decreasing distance than does 
CO(3-2)/CO(2-l), because the CO(3-2) line emission 
is the next to become resolved, as this line is formed more 
inwards in the circumstellar shell, but still further out than 
the higher transitions. 

The stochastic or turbulent velocity determines the in- 
teraction length along the line-of-sight, i.e. the region over 
which the line is formed (see also Sect. I3.ll and l3.2|l . The 
effect of a larger turbulent velocity is different for optically 
thick and optically thin lines. In the optically thin case, 
a change in profile strength may result from changes in 
the line source function in the (near and far) parts of the 
line interaction region, that is added relative to the de- 
fault case. In the optically thick case the relevant source 
function is the one at the location where r « 1, which 
shifts towards the observer when u sto is increased. It may 
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Table 3. Parameters for the temperature profiles with a 
change in slope. See text for details. 
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Fig. 7. Mosaic of diagrams representing the 
CO(4-3)/CO(2-l) ratio versus the CO(3-2)/CO(2-l) 
ratio. Again, the standard model parameters (Table (2J 
were used. Only the temperature profile was varied. From 
the panel in the upper left corner turning clockwise the 
adopted temperature profiles are: i) T(r) oc r~ Q , with a 
indicated in the plot, ii) T(r) is described with a function 
consisting of two power-laws. See text for description, 
iii) T[r) oc 7'~ Q , with negative values of a, indicated 
in the plot, iv) A constant temperature throughout 
the circumstellar shell, where the adopted values are 
indicated in the plot. Again, the predicted line ratios 
are marked with asterisks and the observed line ratios 
with diamonds. The diamonds are not labelled to avoid 
a crowded plot, but can easily be identified using Fig.0] 
The dashed lines indicate equal line ratios. 



therefore differ from the default case. As these effects tend 
in the same direction for all lines (except possibly when 
lines change from optically thin to optically thick) , the line 
ratios are found not to change dramatically when varying 
the turbulent velocity. 

Changing the outer radius has a stronger effect on the 
line ratios, as can be seen in Fig. |SJ When the outer ra- 
dius is increased, more relatively cold gas will be present. 
In this gas mostly the lower rotational levels are popu- 
lated, thus increasing predominantly the CO (2—1) tran- 
sition. The higher the transition, the less it is affected by 
the outer radius. 

The last parameter shown in Fig. H3 is the dust-to-gas 
ratio. The most important effect of adding dust to the 
circumstellar shell is in the population of the rotational 
levels. Continuum emission at 4.6 /zm can be absorbed 
by CO molecules, exciting them from the ground to the 
first vibrational level. They will return to the vibrational 
ground state by spontaneous emission, but preferentially 
to a higher rotational level than they started from. This 
has a non-LTE effect on the level populations, leading to 
variations in both the line strengths and the line ratios. 

We also investigated the effect of the temperature dis- 
tribution. The results are shown in Fig. Various tern- 
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perature profiles have been used in the different panels 
of this figure. The simplest approach is to consider a 
power law T(r) = To(r/R* ) _Q , where a is usually pos- 
itive and has a value ar ound 0.5-0.6 for realistic profiles 
l|justtanont et alJll994 - 

We considered a much broader 
range of a, including negative values and also a constant 
temperature, i.e. a = 0. These cases are of course not a 
true physical representation of the dust shell, but are just 
considered to study the effect of extreme conditions. In 
most cases, we used a temperature at the inner edge of 
To = 2000 K. However, when the power law is shallow 
(small a), the resulting temperature at the outer radius 
would be higher than 25 K if we use the same value for 
To. In that case, we adjusted T such that T(R out ) = 25 K. 
This outer boundary temperature is in the regime of ex- 
citation temperatures of the lower rotational transitions. 
For the negative values of a, the temperature T was as- 
sumed to be 25 K. A number of models with a — has also 
been computed, see the lower left panel of Fig. Adopted 
temperatures are 10, 25, 50 and 100 K. The models with a 
constant temperature or an outwards increasing tempera- 
ture are unrealistic, but we included them in our param- 
eter study, to see if it is possible at all to change the line 
ratios significantly by changing the run of the tempera- 
ture. 

To add to the realism of the models, we composed a 
number of temperature profiles consisting of two power 
laws with different values for a. These profile are inspire d 
by heat balance calculations of lJusttanont et all l)l994f) . 
and are defined as follows: 



T(r) 



T {r/R*Y 



forT> T cx 
for T<T CX 



Five different profiles with a change in slope were 
constructed, where the excitation temperatures of the 
CO(4-3) (T cx = 55.2 K), (3-2) (33.1 K) and (2-1) (16.6 
K) were used to define the position of the change in the 
slope. In all cases a ou t was chosen to be 0.7. For the 
other parameters, the reader is referred to Table [21 The 
resulting line ratios are plotted in the upper right panel of 
Fig. All models with a power law with a slope change 
cluster remarkably close to the CO(4— 3)/CO(2— 1) and 
CO(3-2)/CO(2-l) ratios observed in CRL 2199 and WX 
Psc. The only outlier is profile 2 (see Table [3Jl. Although 
from the various panels in Fig. it seems to be possible 
to explain the observed ratios of WX Psc and CRL 2199, 
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Table 4. Physical parameters of WX Psc. The terminal 
velocity Voo and system velocity of the object wlsr are 
derived from our observations (Table IA.3fl . The other 
parameters are extracted from the lite r ature, the refer- 
ences are: ^van Langevelde et all lll990|) . faofmann et"aT 
" 20011) . ^Justtanont et all lll994l). jZubko fc Elitzui 
2000t) . 5 simb ad . iLe Sidaner fc Le Bertrel l|l996|) . 



Lancon fc Wood! 000' 



parameter value 



ref. 



distance 

VLSR 
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sp. type 
L, 

M» 



0.74 kpc 
20 km s" 1 
+9 km s _1 
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> 5 M Q 
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it is not possible to explain the line ratios of other stars 
of our sample, not even for extreme temperature profiles. 

4.3. A representative case: WX Psc 

In order to investigate the possibilities to explain the inte- 
grated intensities of the CO rotational transitions, we will 
focus on WX Psc. All transitions are observed and de- 
tected. The signal-to-noise ratio is reasonable for all tran- 
sition, except for the CO (6— 5) line. The previous section 
has shown that the line ratios of the lower rotational tran- 
sitions can be explained using power law temperature pro- 
files. In this section we expand our investigations to the 
higher rotational transitions. The observed values for the 
integrated intensities / of the CO (6— 5) and (7—6) tran- 
sition are much lower than the expected values based on 
the standard model described in Sect.0J 

WX Psc is a well studied AGB star with an intermedi- 
at e mass-loss rate. From recent studies, notably the work 
of Irlofmann et all l|200ll) . we have retrieved the physical 
characteristics of the star and the circumstellar environ- 
ment (see Table 0J). These values were used as input pa- 
rameters for our model calculations. For required param- 
eters which are not accurately known, we maintained the 
values of our standard model (Table EJ. We used a stellar 
radius of 5.4 • 10 13 cm, implying a luminosity of 1.3 • 10 
L Q for T cS = 2250 K. 

Determination of the mass-loss rate from the inte- 
grated intensities of the rotational transitions gives an idea 
of the mass- loss history of the AGB star. Fig.[H]shows how 
the integrated line intensities depend on it. In Table the 
mass-loss rates derived from each observed transition are 
listed, while all other parameters were kept fixed. In addi- 
tion, t he gas mass-loss rate , derived from the L — [12 /im] 
colour ^Kemper et al.ll2002|) is given, where a dust-to-gas 
ratio of 1% is assumed. The dust spectral energy distribu- 
tion covers a temperature range of ~ 200 — 800 K, which 
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Fig. 8. Using the known parameters of WX Psc (Table^J 
predicted integrated intensities are given for a large range 
of mass-loss rates. Integrated intensities are plotted in a 
logarithmic scale on the vertical axis, and mass-loss rates 
(Mq yr _1 ) on the horizontal axis, also in logarithmic scale. 
The integrated intensities have been calculated for all lines 
observable with the JCMT and these calculated models 
are indicated with symbols (see legend) . The models have 
been connected with a line. Using the observed integrated 
intensity for a certain line, the mass-loss rate of WX Psc 
can be estimated from this plot (see Table . 



corresponds to a region even more inwards than the CO 
line emission. 

We conclude that constant mass-loss rate models can- 
not explain all of the observed line intensities. Rather, 
it seems that the mass-loss rate varies with the J-level 
under consideration. Specifically, the mass-loss rate cor- 
responding to the CO(2— 1) emission is almost compa- 
rable in strength to the mass-loss rate derived from the 
dust emission (Table |SJ). For the higher rotational transi- 
tions, the derived mass-loss rates go down with increasing 
line frequency, although it perhaps increases slightly again 
for the CO (7— 6) transition. The mass- loss rates deter- 
mined from the high rotational transitions disagree with 
the mass-loss rate derived from the infrared dust emission. 
A difference of at least an order of magnitude occurs al- 
though the regions that are traced by the high rotational 
transitions and the dust emission are closest in tempera- 
ture, and are therefore spatially close together. In general, 
a decreasing mass-loss rate with increasing rotational en- 
ergy level is observed, which is inconsisten t with predic- 
tions based on the superwind model (e.g. iGroenewegenl 
Il994ri Ijusttanont et alJll996t iDelfosse et alJll997|) . The 
results derived here point towards a mass-loss rate de- 
creasing with time, rather than a stratification consistent 
with the onset of a superwind phase. In the next sections 
we will try to explain this discrepancy. 
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Table 5. Mass-loss rates for WX Psc. The values are derived for each observed transition independently. While all 
other parameters of WX Psc were kept constant, the mass-loss rate was determined by fitting the integrated intensity. 
The third column contains the excitation temperature of the corresponding rotational transition. The fourth and fifth 
column represent the outflow distances, i.e. the distance traveled since the gas left the stellar surface, and the last 
column shows the corresponding travel times. In addition to the CO mass-loss determinations, the mass-loss rate 
derived from the L — [12 /im] colour, assuming a dust-to-gas ratio of 0.01, is given ijKemper et alJl2002|) . 
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Fig. 9. Normalized intensity (I(p)p 3 ) as a function of im- 
pact parameter p. The curves for each rotational transition 
are calculated using the mass-loss rate corresponding to 
that transition (see Table EJ. 

4.4. Possible explanations for the inconsistency 
4.4.1. Mass-loss variations? 

In principle, it should be possible to construct a combi- 
nation of a density and temperature profile, such that the 
observed line intensity ratios can be explained. This is 
not possible for a constant mass-loss rate, as becomes ap- 
parent from Fig. so apparently there must also have 
been variations in M. Using the WX Psc model calcula- 
tions, we can derive an impression of the mass-loss history 
using the mass-loss rates listed in Table |SJ For all these 
values we have calculated the region where the respec- 
tive line originates (Fig. 0. This is not done in terms of 
radial distance to the central star, but as a function of 
impact parameter, in which case contributions in line-of- 
sights due to interactions at various radial distances have 
been integrated. Therefore, the values on the horizontal 
axis can not directly be translated to a radial distance to- 
wards the central star, but present a lower limit to this 
distance. In addition, one has to bear in mind that the 



regions from which the various lines originate are not dis- 
tinct, but largely overlap. Some overlap in Fig. 0is due 
to projection effects along the line-of-sight, but a signif- 
icant fraction is due to real physical overlap of the line- 
formation regions. Although all regions are plotted in one 
figure, they do not arise from the same model but are the 
calculated for the corresponding mass-loss rate for each 
line (see Table |SJ). Therefore, it is possible that the CO 
(4—3) seems to originate from a region that is more dis- 
tant from the central star than the region where the CO 
(3—2) line originates, although their excitation tempera- 
tures would suggest otherwise in an outwards decreasing 
temperature profile. Concluding, the mass-loss rates that 
we have determined are only average values for these line 
formation regions. Nevertheless, estimates of the distances 
from the line forming regions towards the central star can 
be derived for the mass-loss rates traced by the observed 
transitions. Using a stellar radius of 5.4 • 10 13 cm and an 
expansion velocity of 20 km s _1 the time elapsed since 
the ejection of the gas from the stellar surface, traced by 
the various transitions can be calculated. The results are 
listed in Table The cycle can be completed by adding 
the dust mass-loss, mostly originating from the region in- 
wards of the CO(7— 6) transition, and thus ejected less 
than 600 years ago. Note that the dust mass-loss rate is 
transferred into a gas mass-loss rate by assuming a dust- 
to-gas mass ratio of 0.01. Actual deviations to this ratio 
imply a different gas mass-loss rate traced by the L — [12 
fim\ colour. From Table we can determine that the in- 
terval between the two maximum mass-loss rates, traced 
by the CO(2— 1) transition and the L — [12 /xm] colour of 
the dust emission, is of the order of ^1000 years. 

Mass-loss variations on such time scales have in 
fact been observed in other evolved stars. Circumstellar 
series of arc-like structures have been interpreted as 
due to mass-loss modulations, notably for C-rich post- 
AGB stars, where the separati on is a measure for the 
time scale of these variations. iKwok et alJ l)l998|) de- 
rive that the separation between arcs observed around 
IRAS 17150—3224 corresponds to a time scale of 
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240 yr (Z?/kpc) (v eKp /10 km s -1 ). For reasonable num- 
bers for the distance and outflow velocity one can de- 
termine that these arcs may be due to mass-loss varia- 
tions on time scales of 200 -1000 yr. A similar time s cale 
(200-800 yr) is derived bv lMauron fe Huggind <|l999|) for 
IRC+10216. The circumstellar arcs around CRL 2688 
(Egg Nebula) are believed to be ejected at 75-200 yr inter- 
vals, assumin g a distance of 1 kp c and an outflow velocity 
of 20 km s- 1 ijsaha.i et alJ liflQaV IRAS LRS spectroscopy 
has shown that hot dust (T > 500 K) is absent around 
a number of AGB stars. This is interpreted as a drop in 
mass-loss rate which occurred ~100 years ago, consistent 
with the sp acing between the ar cs observed around post- 
AGB stars l|Marengo et al.ll200l[) . Hydrodynamic calcula- 
tions considering the gas and dust as partially or com- 
pletely decoupled outflow components resulted in mass- 
loss variations of an order of magnitude at intervals of 
200-350 year fo r partially and 400 year for co mpletely de- 
couple d fluids l|Kimis et al]l200lh / Moreover, iFong et all 
( 200^J| report on the discovery of multiple shells seen in 
CO (1—0) emission around IRC+10216. These shells are 
found to have intershell time scales of 1300-2900 year. The 
circumstellar arcs and molecular shells observed around 
post- AGB stars and the density enhancements emerging 
from hydrodynamic calculations have similar time scales 
to what we derive here for mass-loss variations in the out- 
flow of WX Psc, indicating that the same phenomenon 
may perhaps play a role here. 

Variations in the mass-loss rate of AGB stars have 
already been studied for a long time. It is generally ac- 
cepted that the AGB phase is terminated by the super- 
wind; a phase i n which the mass-loss rate rapidly increases 
llRenzinilll98lHBaud fc HabinJl983h . However, the mass- 
loss rates inferred from the CO line intensities for WX Psc 
decrease with time and are thus opposite to the classical 
superwind model predictions. The thermal pulses associ- 
ated with He-s hell ignition are also though t to cause mass- 
loss variations ijVassiliadis fc Woodlll993j) . As for the su- 
perwind, the behaviour and time scale of these variations 
do not comply with our model predictions. 

4.4.2. A gradient in the turbulent velocity? 

Besides a complex density-temperature profile due to pe- 
riodic mass-loss variations, there may be another way to 
explain the line intensities of the CO rotational transitions 
observed in WX Psc; a gradient in stochastic velocity w sto - 
The stochastic velocity is an important parameter in the 
line formation process (see Sect.|3J). Fig.^Jillustrates this 
for the profiles of one of the rotational transitions for var- 
ious stochastic velocities. 

Analogous to the determination of the mass-loss his- 
tory, it is possible to estimate the variations in the stochas- 
tic velocities traced by the integrated intensities of the 
sequence of rotational transition observed for WX Psc. 
For this purpose, the mass-loss rate was assumed to be 
constant at a rate of 10 -6 M© yr -1 throughout the cir- 
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Fig. 10. The influence on the line profile of CO(3— 2) due 
to variations of the stochastic velocity. The input param- 
eters of our standard model are used (Table |2J) , only the 
turbulent velocity - which in our model is independent of 
r - is varied, in the range from 0.05 - 2.0 km s . 

Table 6. Stochastic velocities for WX Psc. The values are 
derived for each observed transition independently, while 
all other parameters were kept constant. The stochastic 
velocities were determined by fitting the integrated inten- 
sities. 



tracer 




Wsto 

(km s _1 ) 


CO(7- 


-6) 


3.2 ± 0.4 


CO(6- 


-5) 


8 ± 1 


CO(4- 


-3) 


1.0 ± 0.1 


CO(3- 


-2) 


0.24 ± 0.04 


CO(2- 


-1) 


0.16 ± 0.05 



cumstellar outflow. The temperature profile and the other 
parameters were kept the same as the ones used in the 
mass-loss history analysis. The results are listed in Tabled 
Again, the line formation regions and thus the values de- 
rived here are not independent and should be seen as av- 
erages over the formation regions. The observed line in- 
tensities may be explained by a gradient in the stochastic 
velocity if it is lowest in the outer parts of the outflow, 
traced by the low rotational transitions, and has its max- 
imum in the gas traced by the CO(6— 5) transition. 

One has to bear in mind that we derived the stochas- 
tic velocities for one particular mass-loss rate, namely 
1O _6 M0. As pointed out before, the mass-loss rate has 
a considerable effect on the line strengths as well, how- 
ever, the negative gradient will be maintained for other 
choices of M. To explain the observations w s to has to in- 
crease to an unrealistically high maximum of 8 km s _1 in 
the region of CO (6— 5) formation and then decrease again 
to 0.16 km s -1 at the CO(2— 1) formation zone. 

The stochastic velocity can be considered as a compo- 
sition of thermal and turbulent components, according to 
Eq. J5J. The thermal molecular velocities for CO are given 
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by Utherm = \/2kT/mco: where T is the temperature at 
the line formation region, which is usually of the order 
of the excitation temperature. For the lines observed we 
can determine that Wtherm ranges between 0.01 and 0.03 
km s _1 . It is obvious that only a minor fraction of the 
required total stochastic velocity u s to can be explained by 
thermal motion, and also that the observed gradient is not 
sufficiently reproduced by the thermal component. 

The nature of the remaining turbulent velocity «t U rb is 
unknown, but could in part be induced by stellar pulsa- 
tions. These pulsations cause stochastic velocities of 2-5 
km s _1 in the inner parts of the circumstellar shell re- 
quired to start the dust formation process. However, these 
stochastic velocities will damp quick ly and are practically 
absent beyond 100 i?* ijSimislEoOlj) . Hence, if variations 
in the stochastic velocity are important for the CO line 
intensities, the origin of such variations is presently un- 
clear. 

4.4.3. Other factors 

Other factors that could be important include the outflow 
velocity profile, and the geometry. Our adopted outflow 
velocity profile is very simple (see Eq. J2J), but hydrody- 
namical calculations sho w that it may be more complex 
and also time-dependent l(Simis et al.ll200lh . The effect of 
such complex outflow velocity profiles on the CO line pro- 
files has not yet been studied. Perhaps they could serve as 
a source of turbulence. 

Non-spherical winds, e.g. density enhancements in the 
equatorial region, could also play a role in the observed 
line strengths. The observed line profiles would reflect 
such an axi-symmetric geometry, if it exists. Close exam- 
ination of the observed lines (Figs. IA.lllA.2l")l shows that 
their profiles arc similar for all transitions (per source), 
and one can therefore conclude that the regions where the 
lines originate have almost the same velocity structure. 
Apparently there is no change in geometry for the regions 
traced by the various rotational transitions, e.g. a slowly 
outflowing disk traced by the lower transitions and a fast 
polar outflow traced by the higher transitions. Thus this 
possibility most likely can be ruled out as an explana- 
tion for the discrepancy between the observations and the 
model results. Only in case of VY CMa (Fig. \EM thc 
profiles show significant differences between the lower and 
higher rotational transitions. 

5. Concluding remarks 

5.1. CO rotational transitions as mass- loss indicators 

In this work, we presented submillimeter observations of 
various carbon monoxide rotational transitions (CO(7— 6), 
(6-5), (4-3), (3-2), (2-1)) observed in AGB stars 
and red supergiants in various evolutionary states. We 
have attempted to determine the mass-loss history of 
the programme stars by modelling of the observed tran- 
sitions. For the first time the CO(7— 6) and (6—5) 



transitions were used, in addition to lower transitions. 
In this way the gap between the regions in the out- 
flow traced by the gas and that traced by the dust 
emission was largely closed. Many studies in the past 
have focussed on only one or t wo transitions to deter- 
mine the gas mass-loss rate (e.g. JlCnapp fc Morris! Il985l: 
[ Loup et al.1 Il993t lJusttanont et alJ Il994t iGroenewegenl 
I19941J: lJusttanont et alJll996t iDelfosse et al.lll997j) . The 
extension of the data towards higher rotational transi- 
tions clearly demonstrates that determination of a unique 
gas mass-loss rate from a single CO rotational transi- 
tion is highly unreliable. We found that the observed line 
strengths indicate that the outflow has a more complex 
physical structure than was previously assumed. Not a 
superwind, but periodic mass-loss variations comparable 
to the arc-like structures and rings observed around post- 
AGB stars, may possibly account for the observed line 
strengths. Part of the discrepancy could be due to a gra- 
dient in the stochastic velocity as well. 

Independently, another research group has reached 
the same conclusion during the last year. Initially, 
lOlofsson et alJ ((2002) modelled the mass- loss rates of a 
large sample of irregular and semi-regular M-type vari- 
ables by fitting 2, 3 or in one case 4 CO rotational transi- 
tions by assuming a constant mass-loss rate over the last 
1000 years. They derive rates for their sample stars and do 
not report on problems similar to ours, but their figures 2 
and 11 show that the line strength of the higher transitions 
is overestimat ed when this model i s used. In the same vol- 
ume of A&A, ISchoier et ail l|2002h describe a model that 
is able to use periodic mass-loss variations to calculate the 
rotational transitions of CO in C-rich stars. The develop- 
ment of this model is driven by the discovery of mass-loss 
modulations. However, after thorough analysis, they con- 
clude that mass-loss modulations are not important nor 
necessary to explain the CO rotational line profiles . The 
most recent results of lGonzalez Delgado et ail (|2003|) indi- 
cate otherwise, however. When trying to derive the mass 
loss rate of more evolved Miras (i.e. with higher mass- 
loss ra tes than the semi- regulars) . iGonzalez Delgado et alJ 
(2003j) find that a model assuming a constant mass-loss 
rate underestimates the strength of the low transitions. 
This is in principle the same as our result that the high 
transitions are overestimated. 



5.2. Future work 

The work presented here has revealed a much more com- 
plex picture of AGB stellar ejecta than previously as- 
sumed. Additional research is required, which we plan to 
do in the near future. Of particular importance are the 
following issues: 

— First, more observational data should be obtained, in 
particular of high rotational transitions. Our study is 
the first to include the CO(6— 5) and (7—6) transitions 
in the mass-loss rate determinations of three evolved 
stars. In addition, for one object (NML Cyg) observa- 
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tions up to CO(6— 5) were secured. This is not enough 
to draw firm conclusions on the degree of complex- 
ity of the physical structure in the outflows of AGB 
stars, therefore this sample should be enlarged. It is 
important to pay attention to the completeness: if all 
transitions are observed, variations in the important 
physical parameters can be much better constrained. 
In that respect it is also worthwhile to extend the data 
with observations of 13 CO for the lower transitions, 
which provide additional independent constraints on 
the physical conditions. 
— Second, a more realistic representation of the physical 
conditions in the outflow of AGB stars should be used. 
This includes adding a gradient in turbulence and pe- 
riodic mass-loss variations as we have argued in this 
study. In addition, the velocity la w could also be im - 
proved, e.g. following the results of lSimis et al.l l|200ljk 
Although these adjustments will lead to an increase in 
the number of free parameters, it is likely that we will 
be able to use the line profiles to constrain the model 
parameters. This will certainly help in disentangling 
the physical structure of the outflow. 
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Table A. 2. Details of the observations. For each source 
the observed transitions are listed, together with the inte- 
grated observing time in seconds. The correction factor / 
has been applied to our measurements, derived from stan- 
dard measurements. The last column lists the observing 
dates. 
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CO 


J3-2) 


600 


1 


02-Jul-00 


IRC+50137 


CO 


;2-i) 


1800 


0.96 


06-Dec-00 




CO 


J3-2) 


2400 


1.10 


05-Dec-00 


a Ori 


CO 


;3-2) 


5400 


1 


02-Jul-00 


VY CMa 


CO 


;2-i) 


3600 


1 


22-Mar-01 




CO 


J3-2) 


2400 


1.10 


05-Dec-00 




CO 


;6-5) 


8400 


— 


10-Oct-01 




CO 


7-6) 


5400 


— 


09-Oct-01 


a Sco 


CO 


;3-2) 


1800 


1.10 


04-Jul-00 


V438 Oph 


CO 


;2-i) 


1800 


1.01 


03-Sep-02 


AFGL 5379 


CO 


J3-2) 


1200 


1 


17-Apr-00 


VX Sgr 


CO 


;2-i) 


1860 


1 


22-Mar-01 




CO 


;3-2) 


2400 


1 


18-Apr-00 




CO 


;4-3) 


2400 


1 


04-Jul-00 


CRL 2199 


CO 


;2-i) 


1800 


1 


22-Mar-01 




CO 


J3-2) 


1200 


1 


17-Apr-00 




CO 


;4-3) 


8400 


1 


21-Apr-00 


OH 26.5+0.6 


CO 


;2-i) 


1800 


0.95 


22-Mar-01 




CO 


;3-2) 


1200 


1 


17-Apr-00 




CO 


;4-3) 


8400 


1 


21-Apr-00 


OH 30.1-0.7 


CO 


;3-2) 


1020 


1 


06-Jul-00 


OH 32.8-0.3 


CO 


;3-2) 


1200 


1 


18-Apr-00 


OH 44.8-2.3 


CO 


;3-2) 


2400 


1 


07-Jul-00 


IRC+10420 


CO 


;2-i) 


1800 


0.95 


22-Mar-01 




CO 


J3-2) 


1200 


1 


17-Apr-00 




CO 


;4-3) 


2400 


1 


21-Apr-00 


NML Cyg 


CO 


;2-i) 


1800 


0.90 


22-Mar-01 




CO 


J3-2) 


1200 


1 


17-Apr-00 




CO 


;4-3) 


1200 


1 


21-Apr-00 




CO 


;6-5) 


4800 




10-Oct-01 


fx Cep 


CO 


;2-i) 


1800 


1 


20-Sep-02 




CO 


J3-2) 


1200 


1 


18-Apr-00 


IRAS 21554+6204 


CO 


;2-i) 


1920 


1 


20-Sep-02 




CO 


;3-2) 


3600 


1 


07-Jul-00 


OH 104.9+2.4 


CO 


;2-i) 


2400 


0.90 


22-Mar-01 




CO 


;3-2) 


2400 


1 


18-Apr-00 




CO 


;4-3) 


4800 


1.16 


13-Apr-00 








2160 


1 


21-Apr-00 




CO 


;6-5) 


11400 




10-Oct-01 




CO 


7-6) 


2400 




09-Oct-01 
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Table A. 3. Overview of observed line parameters. For each observed transition, the peak intensity (T m b) is measured, 
together with the r.m.s. values (Col. 4) and corresponding bin sizes (Col. 5). Cols. 6 and 7 represent the velocity of 
the object (i>lsr) an d the terminal velocity (v^). The integrated line intensity (Col. 8) is determined by removing 
interstellar absorption and emission from the profile, and integrating the remaining line profile. The accuracy on these 
values is at least 10% for the CO(2-l), (3-2) and (4-3) transitions and at least 30% for the CO(6-5) and (7-6) 
transitions. In some cases the accuracy is deteriorated due to interstellar contamination. 



object 


transition 


Tmb (K) 


r.m.s. (K) 


bin (MHz) 


Vlsr (km s~ x ) 


Vex, (km s" 1 ) 


I (K km s _i ) 


T Cet 


CO 


'2- 


i) 


0.44 


0.056 


0.3125 


+23.1 ± 0.5 


6.7 ± 0.5 


3.9 ± 0.4 


WX Psc 


CO 


'2- 


i) 


2.35 


0.034 


0.6250 


+9.0 ± 0.5 


20.2 ± 0.5 


66 + 7 




CO 


[3- 


2) 


2.91 


0.043 


1.2500 


+9.2 ± 0.5 


20.3 ± 0.5 


82 + 8 




CO 


;4- 


3) 


1.86 


0.090 


0.9375 


+9.0 ± 0.5 


20.6 ± 0.5 


50 + 5 




CO 


;6- 


5) 


0.45 


0.224 


2.5000 


+8.4 ± 1.5 


17 + 3 


10 + 5 




CO 


[7- 


6) 


0.82 


0.378 


3.1250 


+9.8 ± 1.5 


21+3 


23 ± 11 


OH 127.8+0.0 


CO 


[2- 


1) 


0.28 


0.048 


0.6250 


-56 ± 3 


13 + 2 


5.5 ± 1.6 




CO 


[3- 


2) 


0.68 


0.050 


0.6250 


-55 ± 3 


13 + 2 


12 ± 1 




CO 


;4- 


3) 


0.15 


0.057 


1.5625 


complex 


complex 


23 + 2 


o Cet 


CO 


'2- 


1) 


13.80 


0.170 


0.1562 


+46.3 ± 0.1 


8 ± 1 


64 + 6 




CO 


[3- 


2) 


21.61 


0.134 


0.3125 


+46.4 ± 0.1 


8 + 2 


108 ± 11 


IRC+50137 


CO 


'2- 


1) 


1.37 


0.033 


0.4688 


+2.8 ± 0.5 


19.1 ± 0.5 


37 + 4 




CO 


[3- 


2) 


1.44 


0.047 


0.9375 


+3.2 ± 0.5 


18.5 ± 0.5 


39 + 4 


a Ori 


CO 


[3- 


2) 


complex 


0.043 


0.3125 


+3.4 ± 0.5 


15.7 ± 0.5 


50 + 5 


VY CMa 


CO 


'2- 


1) 


complex 


0.033 


0.6250 


+25 ± 3 


47 + 3 


66 + 7 




CO 


[3- 


2) 


3.00 


0.043 


1.2500 


+25 ± 3 


47 + 3 


173 ± 17 




CO 


[6- 


5) 


4.37 


0.469 


3.7500 


+27 ± 2 


48 + 3 


257 ± 77 




CO 


y- 


6) 


7.41 


0.908 


3.7500 


+29 ± 2 


44 + 3 


433 ± 130 


a Sco 


CO 


[3- 


2) 


— 


0.036 


1.8750 




— 


— 


V438 Oph 


CO 


'2- 


1) 


0.24 


0.058 


0.3125 


• ± 0.5 


4.3 ± 1.0 


1.1 ± 0.2 


AFGL 5379 


CO 


[3- 


2) 


2.76 


0.056 


1.2500 


-22.7 ± 1.5 


24 + 2 


84 + 8 


VX Sgr 


CO 


[2- 


1) 


complex 


0.048 


0.4688 


+6.4 ± 1.0 


25 ± 1 


31 + 3 




CO 


[3- 


2) 


2.37 


0.059 


0.6250 


+6.9 ± 1.0 


26 + 1 


97 ± 10 




CO 


'4- 


3) 


1.22 


0.237 


1.2500 


+6.4 ± 0.5 


22 + 2 


42 + 4 


CRL 2199 


CO 


'2- 


1) 


1.16 


0.033 


0.6250 


+33.7 ± 0.2 


17.6 ± 0.5 


26 + 3 




CO 


[3- 


2) 


1.25 


0.034 


1.2500 


+33.5 ± 0.2 


17.9 ± 0.5 


28 + 3 




CO 


[4- 


3) 


0.98 


0.066 


1.2500 


+33.2 ± 0.5 


18 ± 1 


22 + 2 


OH 26.5+0.6 


CO 


'2- 


1) 


contamin. 


0.049 


0.3125 


contamin. 


contamin. 


9 + 3 




CO 


[3- 


2) 


1.05 


0.046 


0.9375 


+26.9 ± 0.5 


18 ± 1 


23 + 2 




CO 


'4- 


3) 


0.86 


0.051 


0.9375 


+27.7 ± 0.3 


15.9 ± 0.5 


19 + 2 


OH 30.1-0.7 


CO 


[3- 


2) 


contamin. 


0.076 


0.9375 


contamin. 


contamin. 


contamin. 


OH 32.8-0.3 


CO 


[3- 


2) 


contamin. 


0.047 


0.9375 


contamin. 


contamin. 


contamin. 


OH 44.8-2.3 


CO 


[3- 


2) 


0.63 


0.047 


0.9375 


-70.3 ± 0.2 


17.7 ± 0.5 


15 + 2 


IRC+10420 


CO 


'2- 


1) 


1.65 


0.034 


0.4688 


+75 ± 1 


43 + 3 


95 ± 10 




CO 


[3- 


2) 


3.23 


0.075 


0.9375 


+75 ± 1 


45 + 3 


180 ± 18 




CO 


'4- 


3) 


2.84 


0.078 


1.8750 


+76.1 ± 0.5 


42 + 2 


150 ± 15 


NML Cyg 


CO 


'2- 


1) 


complex 


0.039 


0.3125 


-2 + 3 


33 + 3 


99 ± 10 




CO 


[3- 


2) 


complex 


0.104 


0.6250 


-2 + 3 


33 + 3 


210 ± 21 




CO 


'4- 


3) 


complex 


0.106 


1.2500 


-1 + 2 


34 + 2 


133 ± 13 




CO 


[6- 


5) 


complex 


0.232 


2.5000 


+2 + 2 


34 + 2 


111 ± 56 


/x Cep 


CO 


'2- 


1) 


complex 


0.025 


1.2500 


+22 ± 2 


33 + 3 


2.5 ± 0.3 




CO 


[3- 


2) 


complex 


0.066 


0.9375 


+21 ± 3 


35 + 3 


14 + 3 


IRAS 21554+6204 


CO 


'2- 


1) 


0.58 


0.041 


0.4688 


-19 ± 1 


19 + 2 


12 ± 1 




CO 


[3- 


2) 


0.49 


0.036 


0.9375 


-19.0 ± 0.5 


18.1 ± 0.5 


10 ± 1 


OH 104.9+2.4 


CO 


'2- 


1) 


0.21 


0.026 


0.6250 


-26 ± 1 


18.6 ± 0.5 


5.4 ± 0.5 




CO 


[3- 


2) 


0.43 


0.043 


0.9375 


-25 ± 1 


18.3 ± 0.5 


11 ± 1 




CO 


'4- 


3) 


0.11 


0.042 


1.8455 


-26 ± 1 


18 ± 1 


3.6 ± 0.5 




CO 


[6- 


5) 




0.224 


3.7500 










CO 


y- 


6) 




1.011 


3.1250 
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Table A. 4. Observed CO line transitions in semi-regular 
variables and AGB stars, obtained from the literature. 
The listed stars are observed in CO(4-3), CO(3-2 ) 



and CO(2— 1). The references are feej^miaji^^dj 
t ] Kerschbaum fc Olofssonl lll999h . c Kn app et al.1 
Yound Jl99fih . iStanek et alJ Jl995l) . " 



2000), 



1998), 



source 


line 




7(K km s _1 ) 


telescope 


ref. 


RV Boo 


CO(4- 


3) 


8.87 


JCMT 


a 




CO(3- 


2) 


6.81 


JCMT 


a 




CO(2- 


1) 


4.47 


JCMT 


a 


X Her 


CO(4- 


3) 


42.17 


JCMT 


b 




CO(3- 


2) 


25.69 


JCMT 


b 




CO(3- 


2) 


13.3 ± 1.3 


CSO 


c 




CO(2- 


1) 


6.4 ± 0.4 


CSO 


c 


R LMi 


CO(4- 


3) 


7.1 


CSO 


d 




CO(3- 


2) 


9.5 


CSO 


d 




CO(2- 


1) 


2.72 ± 0.38 


CSO 


c 


R Hya 


CO(4- 


3) 


46.1 


CSO 


d 




CO(3- 


2) 


37 


CSO 


(1 




CO(3- 


2) 


22.2 ± 2.2 


CSO 


c 




CO(2- 


1) 


4.9 ± 0.5 


CSO 


c 


S Vir 


CO(4- 


3) 


2.3 


CSO 


d 






9"! 

z ) 


z . o 


PSD 


u 




CO(2- 


1) 


0.64 ± 0.25 


CSO 


c 


S CrB 


CO(4- 


3) 


4.9 


CSO 


(1 




CO(3- 


2) 


9.1 


CSO 


(1 




CO(2- 


1) 


2.53 ± 0.51 


CSO 


c 


RU Her 


CO(4- 


3) 


5.7 


CSO 


d 




CO(3- 


2) 


9.4 


CSO 


d 




CO(2- 


1) 


2.3 ± 0.2 


CSO 


c 


X Cyg 


CO(4- 


3) 


52.7 


CSO 


d 




CO(3- 


2) 


63 


CSO 


d 




CO(3- 


2) 


41.5 


CSO 


e 




CO(2- 


1) 


28.8 ± 0.7 


CSO 


c 



-J 



S 
E-, 



-J 




-50 





v (km/s) 



50 



Fig. A.l. JCMT observations of rotational transitions of 
CO observed in WX Psc. See text for details on data re- 
duction and analysis. 
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Fig. A.3. Idem - OH 127.8+0.0 



Fig. A.5. Idem - IRC +50137 
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Fig. A.8. Idem - V438 Oph 
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Fig. A.ll. Idem - VX Sgr 



Fig. A.13. Idem - OH 30.1-0.7 
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Fig. A.15. Idem - OH 32.8-0.3 



Fig. A.17. Idem - IRC +10420 
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Fig. A.20. Idem - IRAS 21554+6204 
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Fig. A.21. Idem - OH 104.9+2.4 



